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Abstract. - The standard expression of the high-temperature Casimir force between perfect 
conductors is obtained by imposing macroscopic boundary conditions on the electromagnetic 
field at metallic interfaces. This force is twice larger than that computed in microscopic classical 
models allowing for charge fiuctuations inside the conductors. We present a direct computation 
of the force between two quantum plasma slabs in the framework of non relativistic quantum 
electrodynamics including quantum and thermal fluctuations of both matter and field. In the 
semi-classical regime, the asymptotic force at large slab separation is identical to that found in 
the above purely classical models, which is therefore the right result. We conclude that when 
calculating the Casimir force at non-zero temperature, fluctuations inside the conductors can 
not be ignored. 



Casimir showed in 1948 [1] that the zero-point energy of the quantum electromagnetic field 
generates an attractive force between two perfectly conducting metallic plates at distance d 
and zero temperature. In his calculation, the microscopic structure of the conductors is not 
taken into account. The latter are merely treated as macroscopic boundary conditions for the 
electromagnetic fields requiring the vanishing of the tangential electric field. This geometrical 
constraint modifies the field eigenmodes depending on d. The d-dependence of the modified 
zero-point energy is the source of the well known Casimir force 

{Ti denotes Planck's constant, c the speed of light). 

The generalisation of Casimir's calculation to thermalized fields was given some years later 
in [2,3], see [4] for a recent account. When the temperature T is different from zero, one can 
form the dimensionless parameter a — (3TThc/d (the ratio of the thermal wave length of the 
photon to the conductors separation; (3 is the inverse temperature). A large value of a (low 
temperature, short separation) characterizes the quantum regime whereas a small value of a 
(high temperature, large separation) yields a purely classical asymptotic result (independent 
of h and c) 

/— i^ + 0(e-n .->0, .>0 (2) 
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where ({s) is the Riemann zeta function. Each field mode is a thermahzed quantum mechan- 
ical oscillator with frequencies obtained from the previously described macroscopic boundary 
conditions. All fluctuations inside the conductors are ignored. We note that in fact, on purely 
dimensional grounds, a term oc must also be proportional to ksT, the only issue being 
the numerical value of the proportionality constant. This issue is the subject of this letter. 

In recent times, a number of works have adressed the question of the incidence of the 
microscopic charge and field fluctuations inside the conductors on the Casimir force [5] , [6] , [7] . 
The considered models are classical : the conductors are represented by slabs (or surfaces) 
containing mobile charges in thermal equilibrium and interacting through the sole Coulomb 
potential. These models all yield the same universal result for the mean electrostatic force 
between the slabs at fixed temperature and large distance 

Universality means that the asymptotic force does not involve any parameter characterizing 
the material constitution of the conductors: particle charges and masses, densities p and slab 
thicknesses. (^) In [5], the authors study a statistical mechanical system of charges confined to 
a plane at distance d of a macroscopic (non fluctuating) planar conductor. In [6], they show 
that replacing the above macroscopic conductor by fluctuating charges docs not alter the 
result (3). We provide in [7] a general derivation of (3) showing that universality is guarantied 
by perfect screening sum rules [8] . 

If one compares the result (3) with (2), one sees that the extrapolation of Casimir's calcu- 
lation to the classical regime is larger by a factor 2 than that obtained in the classical micro- 
scopic models. The two approaches are based on different premises : (2) was derived from the 
frequency spectrum of the full electromagnetic fleld but treating the metals as macroscopic 
bodies without internal structure. One the contrary, the force in (3) is purely electrostatic 
(longitudinal fleld) and it originates from the particle fluctuations inside the conductors. 

This calls for a more complete model that incorporates the dynamical part of the field 
(transverse field) in addition to the internal degrees of freedom of the conductors. A prelimi- 
nary remark is in order: it is well known that classical matter in thermal criuilibrium always 
decouples from the transverse field because of the Bohr-van Leeuwen theorem [9] . It is there- 
fore necessary to treat the conductors' charges quantum mechanically. The complete model 
is formulated as follows. One considers two parallel slabs A and B of surface L^, thickness 
a and at a distance d apart. They contain non relativistic quantum charges (electrons, ions, 
nuclei) with appropriate statistics. The total charge in each slab is taken equal to zero. The 
slabs are immersed in a quantum electromagnetic field, which is itself enclosed into a larger 
box K with sides of length R, L, a. The Hamiltonian of the total finite volume system 
reads in Gaussian units (^) 

i i<j ' ' ^' i 

The sums run on all particles with position and species index 7^; U™''^(7i,ri) is a steep 
external potential that confines the particles in the slabs . It can eventually be taken infinitely 
steep at walls' position implying Dirichlet boundary conditions for the particle wave functions. 

(^)In microscopic conductor models, there is a new energy parameter e^/p~^/*, the mean potential energy, 
so that universality does not follow from a simple dimensional analysis. 

(^)The Pauli coupling terms between spins and magnetic field are not taken into account here. 
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The electromagnetic field is written in the Coulomb (or transverse) gauge so that the 
vector potential A{r) is divergence free and Hq'^'^ is the Hamiltonian of the free radiation 
field. For it wc impose periodic boundary conditions on the faces of the large box K. Hence 
expanding A{r) in the plane waves modes k = (2^^, ^^^) gives the usual formulae 

Mr) = {^y'T.9ik)^^ialxe-''^- + a,,,e^'^-n (5) 



fe,A 

^o"^ = X! ^'^*= '^k,xO'k,\, Wfe = c|A;| (6) 

fc,A 

In (5), efe(A), A = 1,2, are the polarization vectors and g{k), g{0) = 1, is a form factor that 
takes care of ultra-violet divergences. 

Wc suppose that the matter in the slabs is in thermal equilibrium with the radiation field 
and therefore introduce the finite volume free energy of the full system at temperature T 

^R,L,d = -kBT\nT:te-^" (7) 

where the trace Tr = Tr^at Triad is carried over particles' and field's degrees of freedom. The 
force between the slabs by unit surface is now defined by 

/(rf)= hm \im fn,L{d) with /^_^(d) = (8) 
Adding and substracting the free energy of the free photon field in (7) leads to 

^R,L,d = -fcsTln [-^j - keTlnZl^ (9) 

where Z^^'^ is the partition function of the free photon field in the volume K. Since the last 
term of (9) is independent of d, it does not contribute to the force (8). Therefore 

m = ^bT ^ih. -- m (^^^ j (10) 

In principle f{d) yields the Casimir force taking into account quantum and thermal fluctua- 
tions of both matter and field. 

The main result presented in this letter is that, in the semi-classical regime, the dominant 
term of the large distance behaviour of the force (10) is still given by the universal classical be- 
haviour (3). This regime is obtained when the particle thermal wave lengths = ?i(/3/m-y)^/^ 
are much smaller than the slabs' thickness and separation (A^ <C a, d). 

More precisely, the force is of the form 



f{d) = ^^^+n{f3,h,d), where 7e(/3, n, d) = 0(0 (11) 



namely, the quantum corrections included in the remainder 1Z{P, h, d) only occur at the sub- 
dominant order d~^. 

The formalism adapted to the investigation of the high temperature (or semi-classical) 
regime is the Feynman-Kac-Ito path integral representation of the Gibbs weight. In this 
formalism a quantum point particle of species 7 is represented by a closed Brownian path 
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r + A-y^(s), < s < 1, ^(0) — ^(1) — 0, starting at r and of extension A-^: it can be viewed 
as a charged random wire at r. Thus the ensemble of wires can be treated as a classical- hke 
system with phase space points {ri,^i). The wire shape A^^(s) (the quantum fluctuation) 
plays the role of an internal degree of freedom; see [10], section IV, for more details on this 
formalism. Here, for simplicity, we use Maxwell-Boltzmann statistics for the particles. We 
also treat the field classically on the ground that the spacing between the dimensionless energy 
levels (3huJk of the fc field mode become vanishingly small in the high temperature and large 
distance asymptotics (a <C 1). A complete presentation will be found in [11], [12]. 
The Gibbs weight associated to n wires is 

exp(-/3f^e^,e^^y(n,Ci,r,-,Ci) + «X^ / d^jis) ■ Mrj + K^j{s))] (12) 

\ i<j 0=1 ^ ' Jo J 



where 



VM,r,,^j) = I ^^ |,^ + A,,|.(.)-.,-A,..|,-(.)| 

is the Coulomb potential between two wires and the vector potential part a stochastic line 
integral that represents the flux of the magnetic field across the wire. The vector potential 
is itself a random field distributed by the normalized Gaussian thermal weig 
Then the partial trace ( • • • )rad ^ -z^'^r:a.d{^~^^°'"^ ■ ■ ■ ) over the transverse field degrees of 
freedom in (10) is easily performed 








^ ' rad 



i=l 

In (14) Wm is a double stochastic integral 



1 r fju 

e,,e,,W^{nM,r„i,) = ^^^^^^^^ J 5:^.7;(fe, i)G''^(fe)i.(fe,i) (15) 



where 



is the free field covariance and 5i^{k) the transverse Kronecker symbol. (^) In (15), j{k,i) is 
the Fourier transform of the line current j{x, i) = e^^ d^i{s)S{x — Vi — A^.^j(s)) associated 
to the wire ^. One sees that the transverse part of the field gives rise to an eff'ective pairwise 
magnetic interaction Wm that has (up to a factor) the same form as the classical energy of 
a pair of current wires. Its ratio to the Coulomb energy (13) is of the order of kuT divided 
by the rest mass energy of the particles. It accounts for orbital diamagnetic effects, which 



(^)The product in (14) contains the magnetic self energies of the wires. 
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are small in normal conductors. Performing a small k expansion in the integrand of (15) and 
noting that d^{s) = one sees that the large distance behaviour of is dipolar 

e^,e^jWrn{ri,$i,rj,^j) r-. ——1—— f d$i{si) ■ [ d$j{s2) 

X (e7,A^,^i(si) • v..) (e^.A^^^,(s2) • V.J ^ (17) 

Having now identified the basic effective pair interactions between the random wires, 
namely the Coulomb potential V{i,j) (13) and the magnetic potential Wm{i,j) (15), it is 
possible to proceed exactly as in the treatment of classical charged fluids [13] . One sees that 
V{i,j) differs from the genuine classical electrostatic interaction between two charged wires 

•1 rl 1 



Klec(i, j) = / dsi dS2-. , ^ t t \ ^ /■ / M (18) 

Jo Jo In + A^, (si )-rj- X^^ $j (S2) I 

by the quantum-mechanical "equal-time constraint" imposed by the Feynman-Kac formula. 
It is therefore useful to split V{i,j) = 14iec(iji) + Wc{i,j), where 

1 rl 1 



Jo Jo \ri + K^i{si) -rj - Xj^^j{s2)\ 

is the part of V{i,j) due to intrinsic quantum fluctuations {Wc{i,j) vanishes if h is set equal 
to zero). Its large distance behaviour originates from the term bilinear in and ^2 in the 
multipolar expansion of the Coulomb potential in (19). It is dipolar and formally similar to 
that of two electrical dipoles of sizes eiAi^i and e2A2^2- 

{fi ) ) fj 5 ) 

dsi 1 ds2{S{si - S2)) - I) {ejiXj,^iisi) -S/n) {e^jX^j$j{s2) -"^Tj) ] r (20) 

3 Jo P"* ~ fj\ 



/ 

Jo 



Introducing the diagrammatic representation of the correlation func;tions by Mayer graphs, 
we perform the usual resummations of T4iec-chains to sum the Coulomb divergences. This 
provides a short range screened potential ^eiec{hj)j as in the classical Debye-Hiickel mean- 
field theory. Mayer graphs are reorganized in integrable prototype graphs with bonds 

= -/?e^ie^,-^eiec(«,i) (21) 
F^{i, j) = exp[-/3e^,e^, ($eiec(i, j) + WS, j) + W„(i, j))] - 1 + /3e^,e^, $eiec(i, i) (22) 

with the constraint of excluded convolution rule between F{i, j) bonds, namely chains of F 
bonds arc forbidden to avoid double c;ounting of the original Mayer graphs. 

We now sketch the final steps. To obtain the force, one needs to find the asymptotic 
form of the correlation between a wire in A and a wire in B. Set F{i,j) = Fab (Faa) 
when particle i belongs to slab A and particle j belongs to slab B (A), and likewise for 
F^{i,j). Following the methods of [7], one shows that the bond Fab is responsible for the 
universal term — ^(3)/(87r/3(i^) of (11). Some care has to be exercised with the bond F^ab 
that embodies the effect of field and particle quantum fiuctuations through Wm and Wc- It 
has a dipolar long distance behaviour F^{i,j) ~ —(3ej.ej.{Wc{i,j) + Wm{i,j)) ~ \ri — 
that might contribute to the force. In forming the complete correlation function of the two 
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slab system, the bonds Fab and F^ab have to be dressed at their extremities by appropriate 
internal correlations of the individual slabs in conformity with the diagrammatic rules. Thus, 
the complete expressions that enters in the force formula at large separation are of the form 
Gaa * Fab *Gbb and Haa * F^ ab * Hbb- The formation of the slabs' internal correlations 
Gaa and Haa in these terms is not the same because of the excluded convolution rule that 
applies to Fab but not to F^ab- Working out the explicit expressions, one sees that perfect 
screening sum rules in the system of wires applied to Gaa*Fab*Gbb imply the universality 
of the term in (11), but the term Haa*F^ab*Hbb yields no contribution at order 
because of the same sum rules. 

Even without going through the detailed calculations, it is clear from the asymptotic forms 
(17), (20) that the corrections to the electrostatic result (3) due to the quantum nature of the 
charges and the radiation field are controlled by the thermal wave lengths A-y = U^f^Jm^, thus 
small at high-temperature. Because of the Bohr-van Leeuwen theorem, the free energy (7) 
of the complete model continuously approaches that of the corresponding pure electrostatic 
classical system as the A-^ vanish. The force cannot jump by a factor 2 in this limit. 

One must conclude from this analysis that the discrepancy between (2) and (3) is not 
due to the omission of the transverse part of the electromagnetic interaction in the classical 
Coulombic models of refs. [5-7] but should be attribuated to the very fact that fluctuations 
inside the conductors are ignored in the calculation leading to (2). Hence (3) is the correct 
asymptotic form of the high-temperature Casimir force. In other words, the description of 
conductors by mere macroscopic boundary conditions is physically inappropriate whenever 
the effect of thermal fluctuations on the force are considered. 

One the other hand, recent experiments validate the zero temperature formula (1). In 
[14] the authors find an experimental agreement with the value of Casimir force's strength 
7r^?ic/240 to a 15% precision level. This indicates that fluctuations in conductors are drastically 
reduced as the temperature tends to zero and possibly have no more effect on the force at 
T = 0. A full understanding of the cross over from the high temperature regime (11) to the 
zero temperature case together with the role plaid by matter and field fluctuations in the 
conductors is an open problem. 

Finally we like to comment on the Lifshitz versus Schwinger method to take the metallic 
limit in their theories of forces between macroscopic dielectric bodies. In [15] Lifshitz obtained 
the high-temperature large-distance [a. <C 1) force between two dielectric slabs having a static 
dielectric constant e as 



which is easily seen to reduce to (3) in the perfect conductor limit of electrostatics e — > oo. 
In [16], Schwinger et al. have proposed to take the limits in the reverse order, i.e. the perfect 
conductor limit is taken first and the high-temperature large-distance asymptotics afterwards, 
resulting in the value (2). In the light of the preceding considerations, the Lifshitz procedure 
is the right one to recover the high temperature regime for conductors. 




A = 



e+1 



e- 1 



(23) 



* * * 



We thank B. Jancovici for very fruitful discussions on the subject of this work. 



REFERENCES 



[1] Casimir, H. B. C, Proc. K. Ned. Akad. Wet, 51 (1948) 793 



p. R. BuENZLi and Ph. A. Martin: The Casimir force at high temperature 



7 



[2] FiERZ, M., Helv. Phys. Acta, 33 (1960) 855 
[3] Meiira, J., Physica, 37 (1967) 145 

[4] Balian, R. and Duplantier, B., Poincare Seminar 2002 : Vacuum energy-renormalization, 

edited by Duplantier, B. and Rivasseau, R. (Birkhauser Verlag, Basel) 2003 
[5] Forrester, P. J., Jancovici, B. and Tellez, G., J. Stat. Phys., 84 (1996) 359 
[6] Jancovici, B. and Tellez, G., J. Phys. A : Math. Gen., 29 (1996) 1155 
[7] BuENZLi, P. R. and Martin, Ph. A., J. Stat. Phys., 119 (2005) 273 
[8] Martin, Ph. A., Rev. Mod. Phys., 60 (1985) 1075 
[9] Alastuey, a. and Appel, W., Physica A, 276 (2000) 508 
[10] Brydges, D.C. and Martin, Pii. A., J. Stat. Phys., 96 (1999) 1163 

[11] EL BouSTANi, S., BuENZLi, P. R. and Martin, Ph. A., Equilibrium correlations in charged 
fluids coupled to the radiation field, in preparation, ITP, Swiss Federal Institute of Technology 
(EPFL). 

[12] BuENZLi, P. R., and Martin, Ph. A., Microscopic theory of the high-temperature Casimir 
effect, in preparation, ITP, Swiss Federal Institute of Technology (EPFL) 

[13] Hansen, J. P. and McDonald, I. R., Theory of simple liquids, 2"'' ed. (Academic Press, 
London) 1986 

[14] Bressi, G., Carugno, G., Onofrio, R. and Ruoso, G., Phys. Rev. Lett, 88 (2002) 041804 
[15] Lifshitz, E.M., J. Exp. Th. Phys. USSR, 29 (1955) 94 

[16] Schwinger, J., DeRaad Jr., L.L. and Milton, K.A., Ann. Phys., 115 (1978) 1 



